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ABSTRACT 

Because of their neutron-richness and association with supernovae, post- 
explosion protoneutron star winds are thought to be a hkely astrophysical site for 
rapid neutron capture nucleosynthesis (the r-process). However, the most recent 
models of spherical neutrino-driven protoneutron star winds do not produce ro- 
bust r-process nucleosynthesis for 'canonical' neutron stars with a gravitational 
mass of 1.4 and coordinate radius of 10 km. These models fail variously; ei- 
ther the flow entropy is too low, the electron fraction is too high, or the dynamical 
expansion timescale is too long. To date, no models have included the effects of 
an ordered dipole magnetic field. We show that a strong magnetic field can trap 
the outflow in the neutrino heating region, thus leading to much higher mat- 
ter entropy. We estimate both the trapping timescale and the resulting entropy 
amplification. For sufficiently large energy deposition rates, the trapped matter 
emerges dynamically from the region of closed magnetic field lines and escapes to 
infinity. We find that ordered dipoles with surface fields of > 6 x 10^^ G increase 
the asymptotic entropy sufficiently for robust r-process nucleosynthesis. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: mag- 
netic fields — stars: winds, outflows — stars: neutron — supernovae: general 



1. Introduction 

In the r-process, seed nuclei capture neutrons on timescales much shorter than those 
for /5-decay. For sufficiently high neutron-to-seed ratio (> 100), nucleosynthesis proceeds 
to the heaviest nuclei, forming characteristic abundance peaks at A ~80, 130, and 195 
(Burbidge et al. 1957; Cameron 1957). Observational studies of the relative abundance 
pattern of r-process elements in ultra-metal-poor ([Fe/H]< —2.5) halo stars reveal remark- 
able concordance with solar r-process abundances for A > 135 (e.g. Sneden et al. 1996; 
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Burris et al. 2000). These observations suggest an astrophysical r-process site that consis- 
tently produces the same relative abundance in this mass range, at very early times in the 
chemical enrichment history of the galaxy. Although the governing nuclear physics is fairly 
well understood, the astrophysical site for the r-process has yet to be established. Because 
of their intrinsic neutron-richness, the post-explosion outflows expected from protoneutron 
stars (PNSs) in the seconds after successful supernovae (SNe) are thought to be a likely 
candidate site for the r-process (Woosley & Hoffman 1992; Woosley et al. 1994; Burrows, 
Hayes, & Pryxell 1995). 

The characteristics of any potential astrophysical r-process environment that determine 
the resulting nucleosynthesis are the asymptotic entropy (sa; units of /cb baryon"^ through- 
out), the electron fraction (1^^), and the dynamical timescale (rdyn, the density e-folding time 
at T = 0.5 MeV) near the start of the r-process. In general, higher s^, lower Y^, and shorter 

Tdyn, lead to larger ncutron-to-seed ratios and higher maximum A (e.g. Meyer & Brown 
1997). In addition to attaining the requisite Sa, Y^, and Tdyn, a potential r-process site must 
produce a sufficient quantity of r-process ejccta (Mgj) per event, so that the total galactic 
r-process budget is accounted for. For example, if all supernovae produce a robust r-process 
abundance pattern, then the total mass of r-process material that must be produced and 
ejected per SN is M^j ~ 10"^ - IO-^Mq (Qian 2000). 

Recent PNS wind models (not bubble models, as in Woosley et al. 1994) fail to produce 
robust r-process nucleosynthesis up to and beyond the 3"^^ r-process peak for 'canonical' 
neutron stars with M — 1.4 Mq and R — 10 km (Takahashi, Witti, & Janka 1994; Qian & 
Woosley 1996 (QW); Otsuki et al. 2000; Wanajo et al. 2001; Thompson, Burrows, & Meyer 
2001 (TBM); however, see Terasawa et al. 2002). Indeed, given the Y^ and Tdyn derived in 
these models, most efforts fall short of the required Sg, by a factor of 3 — 5. 

To date, protoneutron star wind models have neglected the effects of an ordered dipole 
magnetic field. Nagataki & Kohri (2001) considered the effects of rotation and a monopole- 
like magnetic field on PNS winds in the spirit of Weber & Davis (1967), but were unable 
to explore field strengths above ~10^^ G because of the complicated critical point topologies 
encountered in that steady-state solution. Although QW speculated on the role of magnetic 
fields in their wind solutions, the effects discussed were not quantified. 

Motivated by mult i- dimensional studies of solar coronal hole and helmet streamers, as 
well as by the very high inferred magnetic field strengths of some neutron stars {magnetars, 
Kouveliotou et al. 1999; Duncan & Thompson 1992), we here consider the effects of an 
ordered dipole magnetic field on PNS winds and the composition of their ejecta. In this 
preliminary investigation, we show that strong dipole fields can trap and heat matter in 
regions of closed magnetic field lines and that this matter may be ejected to infinity with 
significantly amplified entropy. Our estimates indicate that for magnetar-strength initial 
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dipole fields, the entropy increase is sufficient for a robust r-process. 

2. Magnetic Protoneutron Star Winds 

Tfie comparison germane for assessing tfie importance of a dipole magnetic field on the 
wind evolution is the ratio f3 — P/{B'^/87t), between the thermal pressure and the magnetic 
energy density. In Fig. 1 we compare P{r) for two wind models with R = 10 km, M = 1.4 M0, 
and Ll°^ = 3.7 x lO^^ergs"^ (labeled, 'High L^') and = 1.8 x lO^^ergs"^ ('Low L^') 
with (dashed lines) for various surface field strengths (Bq). The spherical wind 

profiles shown here are obtained by solving the equations of steady-state general-relativistic 
hydrodynamics in a Schwarzschild metric with a general neutrino energy deposition function 
using a two-point relaxation algorithm on an adaptive radial mesh (as formulated in TBM). 
In Fig. 1, we take B{r) — Bo{RB/r)^, where Rb is a reference radius for the magnetic field 
footpoints. Because of the exponential near- hydrostatic atmosphere, we set Rb — 11km. 
Fig. 1 shows that as the PNS cools and L,y decays, the dynamics of PNS winds can be 
dominated by magnctar-strength magnetic fields. For Bq < 10^"^ G, the magnetic field will 
only dominate the fiow at very late times, when both L,y and the mass outfiow rate (M) 
arc very small. However, for Bq = 10^^ G, the magnetic field dominates even the High L^, 
wind solution for r < i?^ ~ 30 km, where the relationship j3{Rp) = 1 defines i?^. For small 
Bq, such that R13 does not exist, the field is blown into a roughly monopolar configuration 
and the mass flow proceeds to inflnity unencumbered. However, for high Bq, we expect 
the global conflguration to be much like that derived in multi-dimensional studies of solar 
coronal holes and helmet streamers (e.g. Pneumann & Kopp 1971; Steinolfson, Suess, & Wu 
1982; Usmanov et al. 2000; Lionello et al. 2002): open magnetic field lines at the poles and 
a region of closed magnetic field lines (a dead zone, in the spirit of Mestel 1968) at latitudes 
near the magnetic equator. 

In the closed zone the matter is trapped. In the absence of net heating, this configuration 
might be long-lived. In the context of PNS winds, however, we expect the pressure of the 
trapped matter to increase as a result of neutrino energy deposition (g). If P approaches 
B^/Sn in the closed zone, we expect the matter previously trapped to escape dynamically. 
Thus, given /3 < 1 at r (< Rj3) and a local energy deposition rate qp [ergcm~^s~^], the 
amount of time required to increase P to B'^/Sn is 

ru,p [BySn - P]/[qp]. (1) 

Hence, the matter in the closed zone is not trapped permanently, but escapes in Ttrap because 
of neutrino heating. The increase in P of the trapped gas must be attended by an increase 
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in entropy on account of q. If q and T do not change appreciably in Ttrap, Sa increases by 

As ~ grtrap/T. (2) 

If the magnetic field is very weak then Ttrap is very short and As ~ 0. At the other extreme, 
if /3 ^ 1, some of the matter in the closed zone may be permanently trapped (barring MHD 
instabilities, see §4) because as T gets large, neutrino cooling begins to compete with heating 
and q ^ locally. The neutrino cooling rates are stiff functions of T (oc and oc for 
e~p nUg and e~^e~ uu, respectively) and for this reason, T, and, hence, P, cannot 
increase arbitrarily throughout the closed zone. This sets a maximum P increase attainable 
at any radius. For example. Fig. 2 shows T, s, and p (thin solid lines) for the High 
model. Also shown are the maximum temperature achievable at any radius (T*, dotted line) 
for fixed p and the corresponding maximum entropy (s*, dotted line). The thick solid line 
shows [s + As] as computed from eq. (2), for Bq — 5 x 10^^ G. Because [s + As] is greater 
than s* for r < Rq, the material inside Rq is permanently trapped. This defines Rg. Note 
that in all models Rg < Rp. Figure 2 shows that although for very high Bq some material 
may be permanently trapped, the material exterior to Rg and interior to Rp may escape with 
significantly enhanced entropy. In this case, near Rg, Sa is increased from ~ 80 to ~ 600. 



3. Results 

Figure 3 shows Sa versus Tdyn for a 10 km 1.4 PNS, for many different Ly. Each point 
on the thick solid line is a separate spherical steady-state wind model with no magnetic field 
effects. As Li, decreases, the curve evolves from low entropy (~80) and short Tdyn (~4ms) 
to higher entropy (~175) and considerably longer Tdyn (~170ms). One can think of this fine 
as the coohng curve of a PNS composed of snapshots. The dashed fine shows the analytical 
results from Hoffman, Woosley, & Qian (1997) for the Sa required for 3'''^-peak r-process 
nucleosynthesis at fixed and Tdyn- Above this line, for Y^^ = 0.48, 3'''^-peak r-process 
elements can be synthesized. Below it, the neutron-to-sced ratio is too small for the nuclear 
fiow to proceed to A ~ 195. One can see clearly that at all Tdyn, Sa is too low to produce a 
robust r-process for the non-magnetic models. 

The thin solid lines in Fig. 3 are obtained from the non-magnetic models (thick solid 
line) by applying the procedure for calculating Ttrap and As described in §2 at each L^. Lines 
for surface magnetic field strengths {Bq) of 2, 4, 6, and 8 x 10^^ G as well as 1 and 2 x 10^^ G 
are shown. We assume that Tdyn is the same between magnetic and non-magnetic models. 
The increase in enthalpy due to trapping and consideration of the Bernoulli integral indicate 
that this approximation should hold to no better than a factor of 2 — 3. For a given Bq, 
we find that if Rp exists then Rp oc and Ttrap oc i^"^'^, to rough approximation. One 
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can see that for Bq — 2 x 10^ G Rj3 does not exist until Tdyn — 0.075 s. Figure 3 shows that 
the entropy of the flow is profoundly effected by including a strong dipole magnetic field. 
As long as < 0.5, with a surface field <; 6 x 10^^ G, Sa inay be increased sufficiently for 
robust, 3'"'^-peak r-process nucleosynthesis. 

Because the solution to the wind equations constitutes an eigenvalue problem for the 
mass outflow rate (M), we can estimate the amount of r-process material ejected (M^j) by the 
magnetic models by assuming oc t^". For a given Bq, M^- is then the time integral of M 
above the dashed line in Fig. 3. Taking a = 0.9, the entire range of L^, in Fig. 3 corresponds 
to ~ 15s. Given this U{t), MI^{Bq = 2 x 10^^ G)/M^j(5o = 6 x 10^^ G) ~ 50. That is, 
for higher Bq, more r-process material is ejected. Absolute numbers for MJj are misleading 
because M oc R^/^ (QW) and we have not included the PNS radial contraction during the 
early cooling epoch. However, as a first approximation, taking it! = 30 km initially, with 
R oc t~^/^ (until R reaches 10km), with L^, oc t~°'^, and conserving magnetic flux so that 
when R = 10 km, Bq = 5 x 10^^ G, we find that M^- ~ 10^'^ Mq. This number is sensitive to 
Rq, Rj3, and the extent in latitude of the closed zone at the PNS surface, but we conclude 
that although the birth rate of neutron stars with Bq ~ 10^^ G may be small compared to 
the total SN rate, M^j may be large enough to account for the lower rate and accord with 
the total galactic r-process budget (see §1). That MJj is not ejected isotropically may be 
important for the statistics of r-process enrichment in the early galaxy. 



4. Some Open Questions 

Trapping Timescale: Equations (1) and (2) are only applicable if Ttrap is much shorter 
than all other timescales in the problem. We expect global wind quantities to change on the 
timescale for neutrino diffusion and luminosity decay. For example, the Bq — 2 x 10^^ G line 
in Fig. 3 has Ttrap — 0.05s and As ~ 90 for the highest L^, shown (shortest Tdyn)- Because 
q and M do not change significantly over Ttrap, we expect the result presented in Fig. 3 to 
be relatively robust. However, taking the Bq = 6 x 10^"^ G line in Fig. 3, for the lowest Lj, 
(longest Tdyn), we find that Ttrap — H s and As ~ 650. In this case Ttrap is nearly as long as 
the cooling/wind epoch itself and the approximations of tliis paper break down. 

Complex Field Topologies: The magnetic field may be a complex of higher order 
multipoles with a variety of local field strengths in this very early phase of PNS evolution. In 
this case we may consider a distribution of closed loops on the PNS surface. Qualitatively, 
this should not alter the conclusions of this paper - there would simply be a distribution of 
Ttrap and corresponding As locally. For MJj, we need only estimate the covering fraction of 
closed loops. Interestingly, such considerations lead one to consider magnetic reconnection 
as a possible energy deposition mechanism (QW). Naively, including such a term in q may 
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increase the entropy of the affected matter. QW and TBM showed that ~ 10% increases to 
the total energy deposition rate can significantly increase Sg,, but that the radial position of 
an enhancement in q is paramount to determining whether or not the entropy is increased 
asymptotically in the flow. 

MHD Instabilities: The estimate of As in eq. (2) relies on the stability of the closed 
region over a time Ttrap- Given the action of neutrino heating and the geometry of the 
closed zone, it is possible that MHD instabilities might disrupt the system. If the timescale 
for instabilities to grow is much shorter than Ttrap and the instability allows for the rapid 
removal of matter from the closed zone, the entropy enhancements shown in Fig. 3 may be 
significantly decreased. Because of line-tying we expect global interchange instabilities to be 
suppressed (J. Arons, private communication). Ballooning modes at Rp may develop rapidly, 
but allow for only a small amount of matter to escape in Ttrap- Perhaps such instabilities are 
the only way to remove the high entropy matter interior to Rg (see §2). A detailed stability 
analysis, or the solution to the full time-dependent MHD problem is required. 

Rotation & Convection: Similar to MHD instabilities, rapid differential rotation or 
convection may disrupt the global magnetic field configuration on timescales much shorter 
than Ttrap; and thereby compromise the entropy amplifications calculated here. Even the 
shortest significant Ttrap calculated in Fig. 3 is 50 ms in duration, which may be much longer 
than the initial PNS rotation and convective overturn timescales (Duncan & Thompson 1992; 
Thompson & Duncan 1993). 

Spindown: The wind itself may slow the PNS rotation, carrying away angular mo- 
mentum as the PNS contracts and cools. The PNS spins down on a timescale tj = J/ J 
{2,/h)MR^/MRl (e.g. Lamers & CassineUi 1999; Weber & Davis 1967). Taking the 'Low 
L^' solution (M ~ IO-^MqS"^) and Bq = 10^^ G, we find Rp ~ 2000 km and tj ~ 200 s - 
too long in comparison with all other relevant timescales to infiucnce the PNS period signif- 
icantly. For constant i?, we see that for Ly oc and Rp oc L~^-'^^ that tj oct"; the increase 
in Rp as L,^ decreases is insufficient to counter the rapidly decreasing mass fiux. Of course, 
this neglects the coupling between Bq and the rotation of the PNS, which likely exists if the 
magnetic field is generated by either a standard or turbulent/convective dynamo (Thompson 
& Duncan 1993; Thompson & Murray 2001). These estimates for tj, however, assume the 
definition for Rp in §2 and that B oc r~^. Modifications to these assumptions change tj 
significantly. 

5. Summary, Conclusions, & Implications 

For the first time, we consider the effect of a strong dipole magnetic field on protoneutron 
star winds and the resulting nucleosynthesis in their ejecta. We show that magnetar-hke field 
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strengths can dominate the wind dynamics at all times during the PNS cooling epoch (see 
Fig. 1). We argue that at radii close to the PNS, matter may be trapped by closed magnetic 
field lines where (3 < 1. The matter in any closed region is not trapped permanently, but 
escapes on a timescale Ttrap, as a result of neutrino heating and the necessarily commensurate 
increase in both pressure and entropy. As j3 approaches unity in the closed zone, the material 
escapes dynamically. In Fig. 3 we show that for surface magnetic field strengths > 6 x 10^^ G, 
the entropy enhancement engendered by this trapping effect is sufficient to yield robust, 3"^'*- 
peak r-process nucleosynthesis. In this way, the factor of ~ 3 — 5 in entropy previously 
unattained in spherical models of PNS winds may be achieved. That the r-process would 
occur at high entropy (and not, necessarily, at very short Tdyn, as in TBM) in these models 
might help explain how, in the multi-dimensional parameter space of supcrnovae and PNS 
cooling, the same r-process abundance profile above A ~ 135 is observed in both the sim and 
ultra-metal-poor halo stars (see §1). A close look at the models of Meyer & Brown (1997) 
reveals that at s.^ = 500, for Tdyn ~ 0.06, there is considerable leeway in choosing (say 
0.48 to 0.42) so that, approximately, the relative r-process abundances for 135 ^ A < 195 
accord well with those from the sun. For a given Y^, similar flexibility in Tdyn also exists at 
high Sa (also imphed by results of Otsuki et al. 2002). 

Clearly, much more work is required to understand wind emergence, evolution, and 
nucleosynthesis during the cooling epoch of highly magnetized protoneutron stars. This 
preliminary investigation is merely the flrst step in quantifying one of the panoply of possible 
effects. Although the magnitude of the entropy enhancement caused by trapping in closed 
magnetic loops is indicated in Fig. 3, detailed numbers and systematics must await multi- 
dimensional MHD models. That said, however, we here forward the idea that neutrino- 
driven winds from nascent protoneutron stars with magnetar-like surface field strengths are 
responsible for the production of the heavy r-process nuclides we find in nature. 
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the manuscript by the referee, Al Cameron, are also acknowledged. Support for this work 
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for Research in Astronomy, Inc., for NASA, under contract NAS 5-26555. 
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Fig. 1. — Two representative wind pressure profiles (solid lines, logiQ[ergcm^'^]) as a function 
of radius for 'High' and 'Low' neutrino luminosity (L,^) for a neutron star with gravitational 
mass 1.4 M© and coordinate radius R^, — 10 km. Also shown are profiles of magnetic energy 
density (logio[-B^/87r (ergcm~^)], dashed hues) assuming B — Bo{RB/r)^, where the refer- 
ence radius Rb — 11km, for surface magnetic field strengths of 10^^, 10^^, 10^^, and 10^^ G. 

Fig. 2. — T (MeV), s, and p (log^ofg cm~^]), for the 'High L,^' wind solution (thin solid 
lines, compare with Fig. 1). For given p(r), dotted lines show T*{r) and s*(r), the maximum 
temperature and corresponding entropy that can be attained before cooling balances neutrino 
heating (sec §2). For Bq = 5x 10^^ G, the thick solid line shows s + As, where As is computed 
using eq. (2). For this wind model and assumed Bq, Rq ~ 23km and Rjs ~ 125. 

Fig. 3. — Asymptotic entropy (sa) versus dynamical timescale (rdyn) for a collection of PNS 
wind models with R — 10 km and M — 1.4 M© as a function of neutrino luminosity (Li,). 
The thick solid line shows the evolution of the wind solutions in the Sa-Tdyn plane without 
including any magnetic field effects. Thin solid lines show the same curve, but including the 
amplification in Sa due to = 2 x 10^^, 4 x 10^^, 6 x 10^^, 8 x 10^^, 1 x 10^^, and 2 x 10^^ G as 
per eq. (2). The heavy dashed line shows, for = 0.48, the line in the Sa-Tdyn plane above 
which 3'''^-peak r-process nucleosynthesis is likely (from Hoffman, Woosley, & Qian 1997). 
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